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A quadratic extension of REM has been treated. Discussed 
here is the origin of relation of REM to strings and other com- 
plex physical phenomena. Two basic features of the REM 
class of complex phenomena were identified: the double ther- 
modynamic reflection (a hierarchy of free energies) including 
the strong reflection at the upper level (the free energy on the 
order of a logarithm of the degrees of freedom) and the loss 
(complete or partial) of the local symmetry property. Two 
main classes of complex phenomena related to REM are seen: 
the spin glass phase of REM and the boundary the spin glass- 
ferromagnetic phases. Some examples of physics interest are 
analyzed from this viewpoint. 

While quite a number of physical theories (with the 
elements of probability theory) have been successfully 
treated lately, there still remain some unsolved problems, 
the example of which are strings, turbulence, spin glasses, 
2d quantum disorder. The last three ones refer to the 
(loosely defined ) fields of physics of complex systems. 
One can also mention important biophysical systems with 
complex behavior such as the protein folding, the evolu- 
tion. It is noteworthy that however different these 6 fields 
seem, they are associated with REM[1]. It seems reason- 
able, therefore, to establish their relationship to REM 
with a view to better understanding of these theories, 
as well as to comprehend the complex phenomena and, 
if possible, to specify their boundaries, using REM as a 
paradigm of complex systems. 

What is the main feature of REM that makes for its 
power? The human mind is envisaged to be a complex 
system and so the information processing must be an 
essential property of a complex phenomena. From this 
viewpoint REM is unique-it functions as optimal coding 
device by saturating the Shannon limit for highest possi- 
ble rate of error-free coding [2-4] . That is why the model 
may partially resemble even complex adaptive systems of 
the type of biologi cal evolution [5]. 

Recently some interesting data on the relation of REM to 
strings and 2d quantum disorder [6-8] and decaying Burg- 
ers turbulence [9] have been obtained. We are intent to 
clarify further the relation of REM to the quantum Li- 
ouville model on the basis of REM. 

Now consider the standard definition of REM. We have 
M random energy levels distributed according to formula 



P{E)^J — eM-EV'2a) (1) 
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and the spin glass (SG)-paramagnetic (PM) phase transi- 
tion takes place at (3 = w2iS-M point. At complex tem- 



peratures /? = Pi + iP2 there may exist the Lee- Yang- 
Fisher phase (LYF) [10], when f3i < ^. When /? is a 
purely imaginary number, the transition between LYF 
and PM phases occurs at (^2 — Thermodynamically 
REM is equivalent to the p-spin model in the limit of 
large p: 

H = - in..»pSu--s-*p (2) 

l<ji<i2.-<ip<A' 

Here j are the random couplings distributed according 
to normal law. For averaging oilnZ = In^^e"'^'^' in 
distributions Ei, Derrida used the trick: 

/>oo 

lnz = r'(l)+ / lnid[e-*^] (3) 

It is easy to derive this equality for real positive z. Then 
analytically continue that for imaginary z as well. In the 
last case the real part on both the sides of (3) has to be 
taken. To calculate (3), one can factorize the integrand 
in different energies and transform (3) to the form 
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<lnZ>=r'(l)+ / lni(i[/(i)2"] (4) 
where M = 2^, A = , a = N/2 
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/W = 77- / e^-y'^'ny)dy (5) 
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In (5) the integration loop goes round the point from 
the Wright. In this representation the function f{t) can 
be determined for any complex value (t). The main con- 
tribution to (4) in thermodynamic limit ^ oo is given 
by region Ini ^ —oo, so it is possible to consider the sad- 
dle point method in (5). The saddle point is in [—1,0], 
or left wards. When the integration loop is moved to 
encircle the saddle point, different poles of V{y) function 
are intersected in the process, by which means the phases 
of the model are obtained. SG phase is obtained when 
only y — Q pole is is intersected, the Lee- Yang-Fisher 
like phase (LYF) (at complex temperatures) results at 
intersection of y = —2 pole, PM- of y=-l pole. The in- 
tersection of poles to the left oiy — ~2 does not give new 
thermodynamic phases and leads only to singularities in 
finite volume corrections. 

Now consider the normal distribution for Ei with a 
quadratic form of general type (instead of diagonalized 
one in REM), and define the generating function L{t) as: 
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Lit) = ^^=T\ \ - ( dE, 

exp{(-^A,fci?,i?fe)-t^e''^'} (6) 
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In analogy with (3) one can derive 
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\ntdL{t) = -r'(l) + , IT \/ - / 

e«-^E..^^^^^^^)}ln^e^^'5 (7) 



For complex /3 we can assume Re e'^^* to be the def- 
inition of free energy and deduce 
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If case of ReP = , is determined directly by (3) 
and, hence, in (8) instead of L{it) one can consider L{t). 

Introducing the second field instead of complex temper- 
ature, one can construct Z by means of formula: 



^ = ^e/3ii5*cos(/32A) 



(9) 
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where for the fields -D, there is a normal distribution with 
some quadratic form. 

By means of formulas like (8) one can calculate free en- 
ergies Fi 



the free energy in one model is expressed through the 
partition in another by me ans of linear operations. Fi 
models resemble REM. 

When as A the Laplacian is taken, then Li corresponds to 
the quantum Liouville model, L4 to Sine Gordon, -^2,-^3 
are two natural extensions of the Liouville model, de- 
scribing the strings in > 1 dimensions. 
An analogous case of generalized random energy model 
(GREM) assumes rigorous solution. Here a hierarchic 
tree is considered and energy configurations are at the 
end points of the tree. Each point is connected with 
the root of tree via single trajectory. The normally dis- 
tributed random variables arc defined on branches of tree. 
The energy of configuration is determined as a sum of all 
random variables on the trajectory from the root to the 
endpoint corresponding to given configuration. In case 
of GREM the thermodynamic is determined by distribu- 
tion function for one and two levels of energy. So, we 
can try to relate to our Fi the GREM with distribution 
pi(e) =< 6{e - H{s,j)) > - cxp{-^} and /32(e) =< 
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d{e - H{s\j) + H(s^j)) > ^ cxp{--_^_2^- 

for two configurations. It turns out at the consideration 
of p2, that the physics on the whole is determined by the 
character of function 
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Consider the case, when the diagonal elements A^i are 
independent of i. According to Derrida for monotone 
function s' {y) 3 cases are possible. When s'{y) is con- 
stant, then the model is thermodynamically equivalent 
to REM. 

Now let us consider the Liouville model in 2d [6]. For 
L{f) we have 



y D{^{x)exip{~ J d^x(t>{x)'^4>{x)'a-t j d 



} 

(13) 



The coordinate x serves an an index, (^(x) are analogous 
to energy levels Ei and the Laplacian is used for our 
quadratic form. Let us assume the infrared and ultra- 
violet cutoffs as a,L. The number of configurations (of 
different (j){x) is At the calculation for the diagonal 
element of Green function 



(11) 



When the matrix A is large, |^| 00, then Lj and Fj 
should have the same singularities (with respect to rel- 
evant coupling constants). It is not a direct correspon- 
dence, but rather resembles some variant of duality, when 



< 0(f ).^(f ) G{x, x') = ^\n — ^ (14) 

we put |a; — x'l = a. The hierarchy level v stands 
i2-l^Lzl££. It is easy to find, that s{v) = 2v. It is, thus, 
in line with the case of REM according to classification 
of GREM by Derrida [6,7]. In such a case, in the ther- 
modynamic limit one can consider a simple REM with 
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o = G{0, a) and M = ^ (see for mula(l)) and directly 
obtain for the critical point 



(15) 



The correct behavior of s{v) is necessary, but not suf- 
ficient condition for coincidence of the thermodynamics 
of our model (13) in real space with GREM in case of 
s'{v) = const [11]. The authors [6] had highly sound rea- 
sons (in favor of above point)- the coincidence of Green 
functions in both the theories. Though subsequently they 
made an effort to strengthen the arguments, in reality the 
new arguments were no more sound, than that given in 
[6]. 

Now by analogy with (5) we shall try to give another 
representation for L{t) 
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Unfortunately I failed to simplify this expression. 
Let us consider (with A chosen as Laplacian). It cor- 
responds to the to Sine-Gordon model [13]. As it is dual 
to F4, we directly find for the transition point 
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It is connected with PM-LYF transition at /32 = Now 
consider the Liouvillc model at complex temperatures. If 
we analytically continue the formula of DDK 



/3 = A 



s/2b-d- s/T^ 
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into the range d > 1, then at d = 1 we shall have PM-SG 
transition. 

For d > 1 j3 becomes a complex quantity. A possible 
candidate for the analytically continued Liouville model 
is Z/2- In the dual model F2 we find, when the transi- 
tion from SG to LYF phase occurs. As was said before 
(formula (2)), it occurs at /3i = that corresponds to 
d = 19. So besides d = 1, there is another critical di- 
mension for strings at d = 19. In [12] a version of string 
as a candidate for 3d Ising model has been proposed. It 
corresponds to Z/3. 

What is the cause of of the relation of REM to strings? 
Both are representatives of complex systems. 
Let us first set forth the basic principles underlying such 
systems in terms of REM, and then shortly review the 
complex systems. 

It is striking at the phenomenological treatment of com- 
plex systems, that there are two kinds of variables, (spins 



and couplings in case of SG), as well as that at the defini- 
tion the most important models prove to be locally sym- 
metric. Then we shall analyze the properties of REM 
in detail. In SG models one has spins s and couplings 
j. At intermediate stage of consideration of the theory 
of REM, as well as of other models of SG, one comes 
across with the free energy at fixed couplings (the first 
form of free energy) , and then after averaging the second 
form of free energy (the observable one) is encountered. 
Thus, in the beginning Z{j) are calc ulated and then the 
observed free energy for the system as a whole is calcu- 
lated in the form of average < In Z{j) >. In principle, 
instead of In Z one can take the compressed form Z{j) 
as < Z{j)'^ >,!!< 1 [14]. In such a case < Z{jy > cor- 
responds to the partition of the system. In other SG-like 
models (glasses) one can single out the free energies not 
only at the level of the construction of theory, but also 
at the level of observables (macroscopic level). 
REM is distinguished from the other SG models by the 
fact, that the second free energy has a special property, it 
is equal to the logarithm of the number of couplings. The 
effective number in the purest p-spin REM at p = is 
2^, while the free energy (have been scaled to have tran- 
sition at finite temperature) is proportional to N. 
To proceed, let us introduce a concept of reflection as a 
mapping process (with compression in phase space) that 
conserves some essential properties of the initial phase 
space (the possibility to perform macroscopic work-the 
free energy in our case). 

We define the thermodynamic reflection as the one in 
course of which the free energy is determined, and the 
strong reflection, as that where the free energy is the log- 
arithm of the number of particles. 

We define a REM like complex phenomena as the multi- 
ple reflection , the highest hierarchic level of which con- 
tains the strong reflection. 

When at the definition of a model in the framework of 
certain theory we come upon strong reflection and local 
symmetry, we have sound reasons search for its relation 
with REM. Let us define (in somewhat speculative man- 
ner) a concept of local symmetry (or of its loss) for the 
given single sample of SG (instead of an ensemble). 
In case of pure REM nothing observable is changed in 
the PM phase, when the signs of couplings are changed. 
As in case of SG phase it is impossible to identify such 
a symmetry, it is assumed, that the system has lost the 
property of local symmetry. We define the first phase of 
complex systems like this one as REM-SG phase. Let us 
look for other phases of complexity. If in (3) our cou- 
plings in average have ferromagnetic part and there is a 
representation for couplings 



Jo 



(19) 



{jii.Ap ^''^^ normally distributed), in case of sufficiently 
high values of Jo at low temperatures the system is in 
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the ferromagnetic phase with nonzero magnetization < 
Si >= rrii. If our couphng are transformed at fixed 
values of j^, then there takes place a symmetry 

mi^mi'm, ^ jl..i^Vh-Vip (20) 

It is important, that while the symmetry (20) is related 

to the initial one, it is nevertheless a different symme- 
try. When the system moves from the boundary SG-FM 
inside the SG phase, the magnetization m and so the effi- 
ciency of symmetry (20)) decrease. Wc denote this phase 
of complex system as REM- (SG-FM). In the following it 
will be shown, that an analogous situation takes place in 
other complex systems too. In the ordered (simple) phase 
the system has the property of local symmetry (possibly 
in equivalent formulation of theory). This local symme- 
try may be either local gauge symmetry one or higher 
order symmetry (the order of symmetry group is propor- 
tional to the volume of system or is higher by the order 
of magnitude). As the system moves to the higher com- 
plexity phase, it looses the property of local symmetry. 
In principle, one can conceive other case (of weakly com- 
plex phase) resembling PM or LYF phases of REM. A 
bit is similar the transition SG-PM. This corresponds to 
d= 1 case in strings. 

Let us consider now the systems of different complexity 
from our (REM-oriented) viewpoint: presence of logi- 
cally different degrees of freedom; local symmetry (under 
equivalent formulation of problem) ; number of reflections 
and availability of strong refle ctions. 
Turbulence. Is characterized by the presence of Newto- 
nian and dissipative forces, local symmetry (potentiality 
of velocity out of vortices) and multiscaling. It is a good 
candidate for REM- (SG-FM). a short review of diff'erent 
complex systems. 

2d conform theories (c < 1). Here strong single reflec- 
tion and local symmetry is observed. At the calculating 
of correlators REM-PM is encountered. 
Liouville model. Describes the interaction of c < 1 mat- 
ter with 2d gravity. There are strong reflection and local 
symmetry. The model is dual with respect to REM. It 
may be in REM-PM and REM-SG phases. 
Strings. There exist 2d gravity and matter flelds. After 
integration over surfaces with different areas, because of 
the zero mode integration of Laplacian one has double re- 
flections including strong one. REM-SG, REM-PM and 
REM-LYF are possible here. 

Abelian SOC has some local (with very large group) 
symmetries, connected with different ways of toppling. 
Hence, besides the turbulence it is similar to REM-(SG- 
FM) (Though the explicitly local symmetry is not ob- 
served at the formulation of theory). Another probable 
representation of this class is edge of chaos,that is essen- 
tial for the theory of evolution according to S.Kauffman. 
Owing to the presence of multiscaling the availability of 
strong reflection is possible. Though the local symmetry 
here is not observed, REM- (SG-FM) is not excluded. It 



seems proper to remind one of highly interesting intcr- 
pretat ion of this phase (in case of REM) -the system is 
on the threshold of errorless decoding. 
Thus instead of probabilistic (entropic) approach to com- 
plexity I suggest (the free energy approach) to treat the 
complexity in terms of two criteria: 

1. Degree of reflection 

2. Extent of symmetry breaking (loosing). 

If the theory comprises free energy as a logarithmic of 
volume and a local symmetry is present at the dcflnition 
of model, there is a valid chance of relation to REM. 
In REM class of complex phenomena there is an informa- 
tion processing using the lowest level language (any se- 
quence of symbols is allowed, no correlations) due to local 
gauge invariance at the definition of model. In principle 
other classes of complex phenomena with higher com- 
plexity languages may be allowed , as in case of protein 
folding. 
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